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To set up a self-consistent quantum field theory of degenerate systems, the unperturbed state 
should be described by a density matrix instead of a pure state. This increases the combinatorial 
complexity of the many-body equations. Hopf algebraic techniques are used to deal with this 
complexity and show that the Schwinger-Dyson equations are modified in a non-trivial way. The 
hierarchy of Green functions is derived for degenerate systems, and the case of a single electron in 
a two-fold degenerate orbital is calculated in detail. 



The degeneracy of quantum systems is the basis of im- 
portant physical phenomena such as superconductivity, 
magnetism, Bose-Einstein condensation, Jahn- Teller dis- 
tortion, etc. Degenerate systems are also numerous. For 
instance, all systems containing an odd number of elec- 
trons are degenerate by Kramers' theorem 0. Moreover, 
degenerate systems are a subject of active research within 
the density functional approach 0, 0, Q Q ■ 

The quantum field theory of degenerate for quaside- 
generate) systems has a long history 0,0, 0,0 an d devel- 
oped recently into the two-time Green function method 
lCl llll and the covariant-evolution-operator method 
12l UJj |. Still, these theories have a severe drawback: 
either they are not self-consistent or they do not pre- 
serve the symmetry of the system. The antagonism be- 
tween self-consistency and symmetry can be seen on a 
simple example. If the Hamiltonian of a system is spher- 
ically symmetric, the energy eigenstates can be classified 
with angular momentum quantum numbers \LM). How- 
ever, if L > 1, the system is degenerate and no state 
\LM) gives a spherically symmetric charge density [lij . 
So, if that charge density is used to construct a self- 
consistent potential, spherical symmetry is broken. The 
only way to recover the symmetry is to describe the sys- 
tem with the density matrix (2L + l)" 1 YjM \LM){LM\. 
In other words, a self-consistent theory of degenerate sys- 
tems must use a density matrix instead of a pure state. 

In this paper, we set up a quantum field theory of de- 
generate systems based on density matrices. This is a non 
trivial extension of the nondegenerate theory because the 
usual sum over occupied states must be modified, and the 
hierarchy of Green functions acquires additional terms 
that increase considerably the combinatorial complexity 
of the calculations. The problem is solved with two tech- 
nical tools: the nonequilibrium quantum field theory and 
the Hopf algebra of derivations. 

The generating function — In nonequilibrium quantum 
field theory [lfilllfij . the Green functions can be obtained 
by the functional derivative of a generating function Z p 
with respect to double external sources (this is equiva- 
lent to Keldysh' closed-time-path method 0])- In this 
paper we consider nonrelativistic many-body theory and 
the external sources 77+ , ?7_ , 77+ , f)- are fermionic. The 



density matrix is p = ~^2 KL Plk\L){K\ and the generat- 
ing function is 

Z p = tr(pS(fj-,ri-)- 1 S(fj + ,r} + )), (1) 

where S(fj±,V±) = ? 1 exp( - i H int (t)dt + 
2 / V±{ x )i ) { x ) + i> (x)r]±(x)dx) , is the S-matrix, ip(x) is 
the electron field operator and H mt (t) is the interact- 
ing Hamiltonian. According to the functional formula- 
tion of nonequilibrium QFT the generating func- 
tion can be written Z p = c~ lC ' Z®, where D = D + — D- 
and D± is the interaction term H lnt (t)dt where the 
fields are replaced by functional derivatives with respect 
to the external sources fj±,T)±. Finally, Z p is the gen- 
erating function for the noninteracting system: Z p = 
cxp[-iffj{x)G%(x,y)r)(y)dxdy}J2 K LPLKN i j> (L . In this 
equation, rj is a two-dimensional vector with components 
»7_l_ and ?7_, G^{x,y) is the 2x2 matrix 

f -i(0\T(4>(x)^(y))\0) -i(0\^(y)^(x)\0) \ 
\ i(0\i>(x)^(y)\0) -i(0\T*(i>( X W(y))\0) J' 

T* is the anti-time-ordering operator and Nj (L = 
(Jf|:exp [i J fj d (x)ip(x) + ^ (x)r] d (x)dx] :\L), with r\ d = 
r)+ — it- . Notice that N^ L is the generating function for 
the matrix elements of normal products of field opera- 
tors. The proof of these formulas can be found in 
The two technical problems that will be solved in this 
paper are the explicit calculation of N^ L and the manip- 
ulation of e~ lD to recover a hierarchy of Green functions 
for degenerate systems. 

To calculate N^ L , we must define precisely the 
states \K) and \L). The solution of the noninteracting 
Schrodinger equation provides one-electron orbitals u n (r) 
with energy e n . An orbital is called a core orbital if it 
is filled in all states \K) and \L), otherwise, it is called 
a valence orbital. The core orbitals are numbered from 
1 to C, the valence orbitals from 1 to M. There are C 
electrons in the core orbitals and N < M electrons in 
the valence orbitals. For example, in the ion Cr 3+ , there 
are C = 18 core orbitals and N = 3 electrons in the 
M = 10 orbitals of the degenerate 3d shell. The states 
are generated from the vacuum |0) by the action of cre- 
ation operators for core electrons and for valence 
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where i k , j k are valence orbitals (i.e. integers taken in 
the set {1, . . . , M}) ordered so that i% < • • • < ijy and 
ii < • • • < 3n- A lengthy calculation yields 
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with 
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where u n {x) — e~ !E "'u n (r), Un(x) = e le "*-uj l (r) and 
x = (t,r). Notice that a n and a n are anticommuting 
(i.e. Grassmann) variables because the sources are anti- 
commuting variables. 

For the following, we shall need the quantity 

W° p = log Z° p = -i J fj(x)G° Q (x,y) V (y)dxdy 
c 

+ ^2,a k a k + log(p(a,a)), (2) 
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To derive J2J), we used the fact that \og{l+a k a k ) = a k a k , 
because a k and a k are anticommuting variables. For 
the same reason, log(p(a, a)) is a finite polynomial in 
a and a. The term — i f f}{x)GQ{x,y)r){y)dxdy is linear 
in r]± and fj±. The sum over core orbitals is also linear 
in r/± and fj±. So we write log(p(a, a)) = p;(a,a) + 
p c (a,a), where pi(a,a) gathers the terms of log(p(a, a)) 
which are linear in r]± and fj±, and we obtain = 
-i J fj(x)G° p (x,y)r)(y)dxdy + p c (a,a), where G° p (x,y) is 
the sum of all the terms that are linear in rj± and fj+. 

In the standard case of a nondegenerate system Il8j . 
the electron-hole transformation shows that the free 
Green function G°(x,x') is a sum over occupied states 
for t' > t and over empty states for t <t' . However, this 
transformation is usually justified at the first order only, 
by showing that it modifies the free Hamiltonian by a 
constant Equation J2J shows that the electron-hole 
transformation is valid at all orders. In the degenerate 
case, additional terms come from pi(a,a). 

The Hopf algebra of derivations - The hierarchy of 
Green functions will be obtained through a manipula- 
tion of the exponential of functional derivatives e~ tD . 
This manipulation will be made with the help of the 
Hopf superalgebra of derivations. We call derivation 



a functional derivative with respect to rj±(x) or fj±{x). 
We start from the vector space V generated by the basis 
elements S/Sn±(x) and 8/8fj±(x). The exterior Hopf al- 
gebra A(V) is then a standard mathematical object [itl ] 
and we recall here its main properties. The product is 
the composition of functional derivatives (e.g. the prod- 
uct of S/8rj + (x) and 8/8fj-(y) is 8 2 / 8n + (x)Sfj_ (y)). This 
product is anticommutativc. If an element D of A(V) 
is a product of n derivations, we say that its degree is 
deg(Z?) = n and its parity is \D\ = if n is even and 
\D\ — 1 if n is odd. An element D is called even (resp. 
odd) if \D\ = (resp. |D| = 1). The product satisfies the 
identity DD' = (-l)\ D W D ' \ D' D . Thus, an even element 
D commutes with all elements of A(V). We also intro- 
duce a unit 1 such that ID = Dl = D for all D G A(V). 

The main operation that we shall use is the coproduct. 
It is a linear map A from A(V) to A(V) (g> A(V) such that 
Al = 101 and Ad = 9®l + l®9if9isa derivation. 
To define the coproduct of a term D of degree greater 
than 1, we use Sweedler's notation for the coproduct: 
AD = X)-^(i) ® -^(2) an d the recursive definition 

A(DD') = ^(-l)l^)H D a)l(D (1) £>' (1) )®p C2) D; 3) ). 

For example, if d and d' are derivations, 

A(dd') = dd'<Z) l + l®dd' + d <E> d' - d' ® d. 

The convenience of the coproduct stems from the basic 
identity 

D(fg) = J2(- 1 ) ]D " 2)llfl i D M D ^9), (3) 

where / and g are functions of the external sources and 
the fermion fields, deg(/) = n if / is a product of n 
sources or fields and |/| = deg(/) modulo 2. 

The hierarchy of Green functions will be derived from 
an identity involving the reduced coproduct with respect 
to D, denoted by A'. It is defined by A'l = 1 <g> 1, 
A'D = AD - l(g)£)-D(8il and then recursively by 

A'(D n+1 ) = ^(-l^'^WD^D^^D^D^. 

The main property of this reduced coproduct is the fol- 
lowing identity, valid for any even element D of A(V). 
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The hierarchy of Green functions - The one-body 
Green function is G{x,y) = {l/Z p )d'dZ p \ ri ^^, where 
d = i5/8fj±(x) and d' = 5/Sr)±(y). The definition (JTJ 
of Z p , and the fact that tr(p) = 1 yield Z p \ v=fl= o = 1. 
From Z p = e~ w ' Z° p and Z p = e w ° we obtain dZ p = 
e- iD (dZ°) = e- w {8W°)Z°. Noticing that W° is even, 
we operate now with d and obtain 



G(x,y) 



-iD 



((d'dW° p )Z») 
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all quantities being evaluated at r) = 77 = 0. Now we 
apply equations © and and we obtain 

G(x,y) = ^(e-^) (1 , ) (a'aiy p 1 )(e-^) (2 , ) Z p 

-E( e ^)( 1 ')(5W, 1 )(e-^) (2 , ) (9'Z p ), 

where = e^'^W 7 ^. The exponential can be expanded 

G (^) - E H^E(^^)(^>^) 

n=0 

2M-1 

- E ^Ei-i) 10 " 2 '' 1 ^,^ 1 )^)^)^ 5 ) 

n=0 

Equation JSJ is the Schwinger-Dyson equation for degen- 
erate systems. It is the main result of the paper. By fur- 
ther differentiating with respect to external sources and 
using (PJ, we obtain the hierarchy of the n-body Green 
functions. The number of terms of equation J^J is finite, 
but it can be large when M is large. However, the re- 
cursive nature of the coproduct makes it well suited for 
symbolic-algebra packages. Notice that equation (J5J has 
to be used not only for degenerate systems, but more gen- 
erally when the unperturbed state cannot be represented 
by a single Slater determinant. 

The simplest example- As an illustration, we treat in 
detail the simplest case of a single electron in a two-fold 
degenerate level (i.e. N = 1 and M = 2). This exam- 
ple is important because it corresponds to Kramers' de- 
generacy. The Hamiltonian of a nonrelativistic quantum 
system is 

H = [ ^(r)(-A + C/n ( r ))^( r ) dr 
J 2m 

+ L f ^t( r )^)_L^( P '),f,( r )drdr / , (6) 
2 J |r - r'| 

where U n (r) is the nuclear potential. This Hamiltonian 
will be split into two parts as H = Hq + Hi , where Hq is 
the first term of H in JBJ and Hi its second term. Thus 
the interaction is described by the differential operator 
D = D+-D- with 

e 2 r 1 

D± = — I didrdr'- 

2 J |r-r'| 

6* 

6rj± (t, r)5i]± (t, r')Sf}± (t, r')6f}± (t, r) ' 

From this definition we obtain Wp = e~ lD Wp = W?. 

The degenerate levels will be called n = 1 and n = 2. 
The function p(a,a) is p(a,a) = 1 + (l/2)(a 1 ai + 
0:202). Thus log(p(o, a)) = (l/2)(aiai + 0202) — 
( 1/4) Oi 01020:2- The first term is added to the free Green 
function by defining 

G°p{x,y) = G° (x,y)+ipi(x,y) ( 1 ~ 1 ) , 



with pi(x,y) = (l/2)(ui(x)u\(y) + u 2 {x)ul(y)), the sec- 
ond term is p c (a,a) — — (1/4)01010202- Equation (0 is 
now rewritten in terms of Feynman diagrams: 




In this equation, the interacting 1-body Green function 
is G{x, y) = •-«-^-«-• and the same dot with 2n fermion 
lines is the interacting n-body Green function, the free 
Green function is G° p {x,y) — • — - — •, the Coulomb inter- 
action term is >~-« <C , and r(xi,x 2 ,yi,y2) = 

describes the correlation induced by the degeneracy of 
the system: 

r{x 1 ,x 2 ,yi,y2) = - tz , \r- -, — 1—tj — 

or]± {xi)or)± (X2 )orj± {yi)drj± (y 2 ) 

= (-l/4)(ui(xi)u 2 {x 2 ) - u 2 (xi)ui(x 2 )) 
{u\ (yi)4(2/2) ~ ul(yi)u\(y 2 )). 

The first term of Q) corresponds to n = in equation 
(j3J), the next four terms are for n = 1, the next four 
terms for n = 2 and the last term for n = 3. The first 
term and the terms where the initial and final points are 
connected to the white dot (i.e. 3, 4 and 8) come from 
the first line of equation JSJ, the terms 2, 5, 6, 7, 9 and 10 
come from the second line. For a nondegenerate system, 
only the first two terms of the right hand side survive. So 
we see that degeneracy adds a large number of terms to 
the hierarchy. It should be noticed that if equation is 
operated by id/dt — ha (where ft-o = — A/2m + U n ), the 
first term of the right hand side gives S(x — x') 1 the sec- 
ond term remains and all the other terms vanish because 
(id/dt — ho)ui(x) = 0. In other words, the differential 
form of the hierarchy is the same for degenerate and non- 
degenerate systems, but the integral form (|7| is different. 
So the additional terms can be seen as boundary condi- 
tions due to the correlation of the system induced by its 
degeneracy. Finally, it should be stressed that this hier- 
archy is non perturbative. All Green functions involved 



4 



in equation Q are interacting. However, this hierarchy 
can be used to derive a perturbative expansion of the 
Green function (see, e.g. |2l|h 

To solve equation J7J, the hierarchy must be closed 
by using one of the standard approximations such as the 
random-phase approximation, the GW-approximation or 
the Bcthe-Salpeter equation. 

Self-consistency - The previous equations were not 
self-consistent. To obtain self-consistent equations the 
Hamiltonian H will be split into two parts as H = 
Hq + Hu with H =j ^t( r )(_ A. + Un (r) - V»)^(r)dr 
and Hi = H — Hq. The potential V(r) will be determined 
self-consistently. The presence of this external field adds 
the term D± to D±, where 

D v + = / dtdrV(r)- — - — - — ; — -, 

and the differential operator describing the interaction is 
now D SCF = D + +D v + -D--D v _, This new interaction 
does not modify Wp = W9, but it adds the following 
terms to the right-hand side of equation 




The condition of self-consistency is now that the po- 
tential V(r) is generated by the electron charge density 

of the system or x = — • — in terms of Feynman 

diagrams. 

In this paper we showed how to calculate the Green 
function of a degenerate system, and we derived the cor- 
responding hierarchy of Green functions. An advantage 
of our approach is that it preserves the symmetry of the 
Hamiltonian if the closure of the hierarchy is done prop- 
erly. Recall that this is not possible if we use states in- 
stead of density matrices. 

Equation (J5J is written in terms of unconnected Green 
functions. By defining a second type of reduced coprod- 
uct, it is possible to write the hierarchy of connected 



Green functions. In the near future, several develop- 
ments of the present formalism will be presented. Firstly, 
if we add photon fields, the same formalism can be used 
to derive the Schwinger-Dyson equations for the Green 
functions of a degenerate systems in quantum electrody- 
namics; this will be applied to the nonperturbative cal- 
culation of a degenerate atomic system. Secondly, in the 
derivation of equation J£J, no special form was assumed 
for the density matrix p. So the present formalism can 
be used also for out-of-equilibrium quantum field theory. 
In particular, we can calculate the energy of the inter- 
acting system corresponding to a general density matrix 
p. If we minimize the total energy with respect to p, we 
obtain a set of equations that unify the Green function 
and diagonalization methods of many-body theory. 

I gratefully acknowledge the fruitful comments by Mat- 
teo Calandra, Philippe Sainctavit and Francesco Mauri. 
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